Correlation effects in Ni 3d states of LaNiPO 
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The electronic structure of the new superconducting material LaNiPO experimentally probed by 
soft X-ray spectroscopy and theoretically calculated by the combination of local density approxima- 
tion with Dynamical Mean-Field Theory (LDA-I-DMFT) are compared herein. We have measured 
the Ni L2,3 X-ray emission (XES) and absorption (XAS) spectra which probe the occupied and 
unoccupied the Ni 3d states, respectively. In LaNiPO, the Ni 3d states are strongly renormalized by 
dynamical correlations and shifted about 1.5 eV lower in the valence band than the corresponding 
Fe 3d states in LaFeAsO. We further obtain a lower Hubbard band at -9 eV below the Fermi level 
in LaNiPO which bears striking resemblance to the lower Hubbard band in the correlated oxide 
NiO, while no such band is observed in LaFeAsO. These results are also supported by the intensity 
ratio between the transition metal L2 and L3 bands measured experimentally to be higher in La- 
NiPO than in LaFeAsO, indicating the presence of the stronger electron correlations in the Ni 3d 
states in LaNiPO in comparison with the Fe 3d states in LaFeAsO. These findings are in accordance 
with resonantly excited transition metal L3 X-ray emission spectra which probe occupied metal 
3d-states and show the appearance of the lower Hubbard band in LaNiPO and NiO and its absence 
in LaFeAsO. 

PACS numbers: 71.27.-(-a, 74.25. Jb, 78.70.En 



I. INTRODUCTION 

The superconductivity of the quaternary transition 
metal oxyphosphides LaFePO [T^ = 3.2 K) and LaNiPO 
(Tc = 3.0 - 4.3 K) was discovered recentljii"— although 
these compounds have been under study for more than 
ten years."* In spite of the relatively lovif superconducting 
transition temperatures these materials are important 
because they triggered the extensive search for supercon- 
ductivity in oxyarsenides LnMeAsO (where Ln — La, 
Ce, Pr, Nd, Sm, Gd; Me = 3d metals). This search has 
proven to be remarkably successful, reaching TcS of up 
to 55 K^"— and strong upper critical fields Hc2 of up to 
100 T^ 

Despite the focus on oxyarsenides, Ni and Fe oxyphos- 
phides still attract attentioni^"— because while they 
share the same bi-layer structure as LnMeAsO mate- 
rials, the mechanism for superconductivity appears to be 
different. In particular the value of the total electron- 
phonon coupling constant A in LaFeAsO is much lower 
than in conventional electron-phonon coupling supercon- 
ductors, for example, compare the A = 0.21 of LaFeAsOii 
with the A = 0.44 of Al (where = 1.3 K for Al), and 
even the inclusion of multiband effects fails to explain 
the observed Tc of 26 K.** For LaNiPO the coupling 
constant is more then two times higher (A = 0.58) and 
the superconducting properties can be described within 
the Migdal-Eliashberg theoryji^ Herein we interpret soft 
X-ray emission and absorption spectra of LaNiPO and 



LaFeAsOi^ and compare the measurements with our lo- 
cal density approximation with Dynamical Mean Field 
Theory (LDA-I-DMFT) electronic structure calculations 
to investigate the similarities and differences between 
these two types of superconductors. To assist in inves- 
tigating these materials, electronic structure calculations 
of LaFePO^" (T^ = 3.2 Ki) and Ni02i were used. 

LaFeAsO, unlike LaNiPO and LaFePO, is not super- 
conducting unless dopingS^i^ or high pressure2ii2^ is ap- 
plied to suppress the magnetic transition temperature 1^ 
Since these materials share the same basic ambient crys- 
tal structure and atomic constituents, yet exhibit differ- 
ent low-temperature properties, a basic study of the am- 
bient electronic structure of these three materials is of 
interest, especially since the bulk electronic structure of 
these materials is insensitive to temperature or magnetic 
phase changes 1^ 



II. EXPERIMENTAL AND CALCULATION 
DETAILS 

Single crystals of LaNiPO were synthesized by heat- 
ing a mixture of 375.0 mg La (99.9%, Smart Elements), 
201.7 mg NiO (99.99%, Sigma- Aldrich) and 83.6 mg P 
(red, 99%, Sigma- Aldrich) with 2000 mg Sn (99.99%, 
Alfa Aesar) in an alumina crucible, which was sealed 
in a silica tube under an atmosphere of purified argon. 
The sample was heated to 1173 K at a rate of 40 K/h, 



kept at this temperature for 10 days and slowly cooled 
down to room temperature at a rate of 3 K/h. The cru- 
cible was smashed and the tin bar dissolved in 6 M HCl 
at room temperature. The remaining sample consisted 
of single crystals of LaNiPO beside small amounts of 
LaNi2P2 (7%), NisSnP (4%) and Ni3Sn4 (< 1%). Fur- 
ther attempts to optimize the synthesis conditions with 
regard to reaction temperature or duration were unsuc- 
cessful. Samples prepared directly from the starting ma- 
terial without tin flux yielded only small amounts of La- 
NiPO with LaNi2P2 as the main product. For details of 
preparation see Ref. The soft X-ray absorption and 
emission measurements of the metal L2,3 edges were per- 
formed at the soft X-ray fluorescence endstation of Beam- 
line 8.0.1 at the Advanced Light Source in the Lawrence 
Berkeley National Laboratory. ■^^ The endstation uses a 
Rowland circle geometry X-ray spectrometer with spher- 
ical gratings and an area-sensitive multichannel detec- 
tor. We measured the resonant and non-resonant Ni L2,3 
(3d, 4s — >■ 2p transition) X-ray emission spectra (XES) for 
LaNiPO. Additional non-resonant XES measurements of 
the Ni L2,3 edges of Ni metal foil and NiO were obtained 
as reference standard. The instrumental resolving power 
(E/AE) for emission measurements was about 10'^. The 
X-ray absorption spectra (XAS) were measured in total 
electron yield (TEY) mode for the Ni L2,3 edges. The in- 
strumental resolving power (E/AE) for absorption mea- 
surements was about 5 x 10'^. All absorption spectra 
were normalized to the incident photon current using a 
highly transparent gold mesh in front of the sample to 
correct for intensity fluctuations in the incoming photon 
beam. The excitation energies for the Ni ^2,3 resonant 
X-ray emission spectra were determined from the XAS 
spectra and the energies were selected at the L3 and L2 
thresholds. 

Electronic structure calculations were performed 
within the pseudopotential plane-wave method PWSCF, 
as implemented in the Quantum ESPRESSO pack- 
age)^ We used the generalized gradient approxi- 
mation in the Perdew-Burke-Ernzerhof version-- for 
the exchange-correlation potential in the Rappe-Rabe- 
Kaxiras-Joannopoulos form<^ The Brillouin zone inte- 
gration was performed with a 15 x 15 x 15 k-point grid. 
A kinetic-energy cutoff of 45 Ry was employed for the 
plane-wave expansion of the electronic states. The ex- 
perimentally determined lattice parameters and internal 
atom positions of LaNiPO (a = 4.0461 1, c = 8.100 Af- 
were used. 

To include dynamical correlation effects in the 3d 
shell of Ni, we performed the LDA-I-DMET^^ calcula- 
tions for LaNiPO. Follo wing the Wannier function pro- 
jection procedure of Ref. l32l we constructed an effective 
Hloa Hamiltonian and then used it to solve the Dy- 
namical Mean-Field Theory (DMFT)^ self-consistency 
equations. The H^oA Hamiltonian contained 22 bands 
due to five Ni 3d, three O 2p, and three P 3p orbitals per 
formula unit, projected in a single energy window that 
explicitly takes into account the hybridization between p 
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FIG. 1: (Color online) Total and partial densities of states for 
LaNiPO (in comparison with LaFeAsO from Ref. 36) obtained 
within the density functional theory (DFT) calculations. The 
dashed lines in the Ni, Fe 4s, 3d DOS refer to the metal 
4s states magnified by a factor of 10. La does not have any 
significant contribution to the valence band, and is not shown 
here. 



and d electrons)^ 

The DMFT auxiliary impurity problem was solved 
by the hybridization function expansion Continuous- 
Time Quantum Monte-Carlo methodi^ The elements 
of Coulomb interaction matrix were parameterized by 
U and J parameters;^ We used interaction parameters 
U = 8 eV and J = 1 eV for LaNiPO similar to the val- 
ues obtained in Ref. [2l|. Calculations were performed 
in the paramagnetic state at the inverse temperature 
f3 — 1/T — 20 eV^^. The real- axis self-energy needed 
to calculate spectral functions was obtained by the Fade 



approximant 
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III. RESULTS AND DISCUSSION 

The calculated noncorrelated electronic structure of 
LaNiPO is shown in Fig.[T]in comparison with the struc- 
ture of LaFeAsO. "^^ These calculations are in agreement 
with other DOS calculations available to da.teJ^^ In all 
cases the far bottom of the valence band (-11 eV) con- 
sists of P 3s or As 4s. The top of the valence band 
(-2 eV to eV) consists almost solely of metal 3d states 
in both cases. Between -2 and -4 eV there is strong 
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FIG. 2: (Color online) Densities of states for Ni 3d and 
Fe 3d orbitals obtained within DFT (filled areas) and the 
LDA+DMFT total 3d spectral functions (sohd hnes). Data 
for LaFeAsO^ and NiO^i are given for comparison. 
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hybridization between the O 2p and Ni 3d states in La- 
NiPO; in LaFeAsO there are far fewer Fe 3d states in 
this region indicating much weaker hybridization. The 
situation is the same with LaFePO;^ LaNiPO also has a 
reduction in metal 3d states and total states at the Fermi 
level compared to LaFeAsO and LaFePO. This may ex- 
plain why Ni-based superconductors have lower Tc values 
than FeAs-based superconductors. For all compounds 
the P 3s, 3p and As 4s, 4p states, respectively, occupy 
the same basic region in the valence band and do not con- 
tribute significantly to the Fermi level. The La 5p states 
are identical for both compounds, they have atomic-like 
character and do not contribute to the valence band; they 
are not shown here. 

The spectral function from the LDA-I-DMFT calcula- 
tion for the Ni 3c? states of LaNiPO is shown in Fig. [5] 
(middle panel). In the energy interval from -1 to 1 eV 
near the Fermi energy the 3d spectral function is close to 
the noncorrelated LDA density of states. However, below 
these energies the spectral function is substantially renor- 
malized with the formation of a strong peak at -1.5 eV, 
and the appearance of a lower Hubbard band: the broad 
peak centered at -9.1 eV. Thus this picture resembles the 
LDA-fDMFT results for LaFeAsO^'' (upper panel) only 
for the energies above -6 eV, since in LaFeAsO the lower 
Hubbard band was not found. In NiO, however, a simi- 
lar broad peak centered at -10 eV is obtained, as shown 
in the lower panel of Fig. [5J This lower Hubbard band 
is evidence for strong correlations,-^ and this is a clear 
indication of strong correlations in LaNiPO, similar to 
mOM The comparison of the LDA+DMFT calculation 
for Ni 3d states of LaNiPO and NiO and Fe 3d states 



FIG. 3: (Color online) Comparison of the LDA-I-DMFT total 
3d s pec tral functions (LaFeAsO from Ref. and NiO from 
Ref. '21') with the resonantly excited Ni L3 and Fe Ls X-ray 
emission spectra. 



of LaFeAsO with the resonantly excited Ni L3 and Fe 
L3 X-ray emission spectra which probe occupied Me 3d 
states is presented in Fig. [31 The occurrence of the lower 
Hubbard band in LaNiPO and NiO and its absence in 
LaFeAsO is confirmed by experimental XES spectra. 

The soft X-ray metal (Ni, Fe) L2.3 spectra are shown in 
Fig.m The metal £2,3 XES indicate two main bands sep- 
arated by the spin-orbit splitting of the metal 2p states. 
The lower intensity high energy band corresponds to the 
L2 emission line (3d, 4s — 2pi/2 transitions), and the 
higher intensity low energy band corresponds to the L3 
emission line (3d, 4s — >■ 2^3/2 transitions). The resonant 
L2 and L3 XES (curves b and c in the bottom panels, 
respectively) have the same basic shape. The lack of res- 
onant features indicates that the spectra primarily mea- 
suring the partial occupied DOS rather than multiplet or 
inelastic scattering effects. Note that the La M4.5 XES 
appears below the Ni L3 emission line in the resonant 
Ni L3 spectrum. The metal L2,3 XAS are presented in 
the top panels of Fig. [H According to dipole selection 
rules (A/ = ± 1) they correspond to the excitation of 
metal 2p-core level electrons into unoccupied 3d states. 
Unfortunately these spectra can not probe the unoccu- 
pied 3d DOS directly because the core-hole causes an in- 
creased effective nuclear charge distorting the local DOS 
levels. Further, simulating L2.3 XAS requires consider- 
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FIG. 4: Summary of spectra for LaNiPO (left panel) and 
LaFeAsO (right panel). The upper panels show the metal 
L2,3 XAS (in TEY mode), the lower panels the resonant and 
non-resonant metal L2,3 XES. The excitation energies are in- 
dicated by arrows in the XAS plots. 



ing multiplet splitting, hybridization, and crystal field 
effects. One such simulation was recently conducted for 
LaFeAsO in Ref. |3^, to our knowledge no similar simu- 
lation of LaNiPO exists. Therefore we include the metal 
L2,3 XAS only for completeness. Resonantly excited Ni 
L3 XES of LaNiPO (curve c) shows the presence of La 
M4^5 XES because excitation energy in this case is very 
close to resonant excitation of La M-emission spectra. 

The ratio of the integral-intensity of the metal L2 and 
L3 peaks (the 1{L2) /liL^) ratio) for LaFeAsO is roughly 
the same as that of metallic Fe, and quite different from 
that of strongly correlated FeO (see Fig. [5l right side, 
bottom panel). In a free atom, the I(L2)/I(L3) ratio 
should be equal to 1/2 as the ratio is based solely on 
the statistical population of the 2pi/2 and 2p3/2 lev- 
els. In metals the radiationless L2L^Mi^cj Coster-Kronig 
(C-K) transitions greatly reduce the 1{L2) /1{L^) ratio^^ 
and the I(L2)/I(^2) ratio can be used as a measure for 
the electron correlation strength of a transition metal 
compound^ (see Fig. [5l right side, top panel). The 
full width at half maximum (FWHM) of the L3 band in 
LaFeAsO is again closer to that of metallic Fe than FeO 
(see Fig. [Sj right side, middle panel). While this does 
not directly prove anything, it suggests that the Fe 3d 
electronic structure of LaFeAsO may be similar to that 
of metallic Fe. The shape and statistics of the Fe L2,3 
XES indicate that the Fe 3c? states in LaFeAsO are not 
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FIG. 5: Comparison of non-resonant metal L2,s XES of NiO, 
LaNiPO, and Ni (left side panels) and FeO, LaFeAsO, and 
Fe (right side panels) from Ref. [Igj. The \(L2)/l{L'i) ratios 
for each system are shown in the top panels, and the full 
width at half maximum (FWHM) of the L3 bands are shown 
in the middle panels. The metal L2,3 XES are shown in the 
bottom panels, for easy reference. The I(L2)/I(i3) ratios 
were calculated by taking the quotient of the integrals of the 
Z/2 and Z/3 bands. 



strongly correlated. 

In contrast to LaFeAsO, the I(L2)/I(L3) ratio for La- 
NiPO is much greater than that of Ni metal, as is the 
FWHM of the LaNiPO L3 band (see Fig. [3 left side, 
bottom panel). Indeed, the 1{L2)/1{L^) ratio and L3 
FWHM for LaNiPO (see Fig. [3 left side, top and mid- 
dle panels) are rather close to those of correlated NiO. 
Since the transition metal 1{L2) /1{L^) ratio is over 50% 
greater in LaNiPO than LaFeAsO, and since in NiO is 
comparable to FeO in terms of "correlation strength" ^ 
the Ni 3d states are more correlated than the LaFeAsO 
Fe 3d states. 



IV. CONCLUSIONS 

We have studied the electronic structure of LaNiPO ex- 
cited by synchrotron soft X-ray emission and absorption 
spectroscopy and obtained the theoretical spectral func- 
tions within the combination of local density approxima- 
tion with Dynamical Mean-Field Theory (LDA-^DMFT). 
We conclude that the Ni 3d states of LaNiPO reside 
deeper in the valence band than the Fe 3d states of 
LaFeAsO. The greater occupation in the metal 3d bands 
in LaNiPO reduces the density of the states at the Fermi 
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level and increases the hybridization with O 2p states 
compared to those in LaFeAsO. Accounting for dynam- 
ical correlation in the Ni 3d states of LaNiPO results in 
the renormalization of the states below the Fermi energy 
and the formation of the lower Hubbard band centered 
at -9 eV, similar to NiO, but in contrast to LaFeAsO. 
The I(Z/2)/I(i'3) ratio is much higher in LaNiPO than in 
LaFeAsO, indicating the Ni 3d states of LaNiPO have 
stronger electron correlations than the Fe 3d states of 
LaFeAsO. 
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